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ABSTRACT 


Some Nonlinear Feedback Control Systems under high gain con- 
diticss exhibit the phenomenon of subhsrmonic instability, or contain 
subhzrmoics of the fundamental cutput frequency. A general discussion 
of subhsarminics in monlinear systems is foilowed by an investigation 
of ripple instability in a forced limit-cycling voltage regulator con- 
taiming a thyristur or SCk bridge utilizing an ON-OFF switching 
scheme. 


A digital cemputer prvcgram 1s used tu simulate the dynamic 


resp use uf the system under different loading conditions and for 


different reference vultage levels. 
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ee Introduct Lon 


There are many physical applications where a very stable direct 
current power supply is desirable. This is usually accomplished by 
taking 4a multiphase slcternmating current signal, full wave reccifying 
it and then passing it through whst is essentially a low pass filter. 
The output then 1s some DC level cn which rides a small value of 
"ripple". The ripple frequency is some multiple of the source input 
sigaal and is due to the filter's being non-ideal, 

For m°st applications the amplitude of the ripple voltage is 
not so great that the power supply cannot be used and considered 
tu be a stable DC source, 

However, with the advent of the computer and switches that 
need only millivelts of fluctuations to switch, it is desirable to 
have power supplies that are more tightly regulated and effect a 
fast reccvery from fluctuations 1n source supply and noise pertur- 
bsatiuns, This scheme is usually accomplished with higher loup gain. 

When the loop gain increases to higher values, however, it 
has been witnessed thac the amplitude of the ripple increases and 
contains components of frequencies that are lower than those wit- 
nessed under Lower gain conditions, 

Power supplies that employ silicon cuntrolled rectifiers 
(thyristcers) which are of the pha@se-cuontrol type are known to contain 
these subharmonics under high gain Micsticus..- 

The nature uf this paper is to investigate and verify the 
eccurrence cf subharmonics in a forced limit-cycling regulator em- 
plcoyving silicon controlled rectifiers that are not phase fired. 


While verifying the fact chat they do exist, an investigation 1s 


made in an attempt to determine how these subharmonics are caused, 
what frequencies are present, and under what cunditions they occur. 
The ultimate gcal of such investigation, of course, would be 
t> disccver means two predict the cccarrence of gubharmonics, their 
frequercies and amplitudes, with hopes of designing means to com- 
pensste the system so that these subharmenics do not present too 


great of a degrseding effect on the output of che regulated power 


supply. 


2. A general study of subharmeonics in @ nonlinear system, 

If a linear system is forced by a sinusoidal input, the 
resulting steady-state output 1s also a sinusoid, However, a non- 
linesr system driven by a sinusoid dces not produce a sinusoidal 
output, but the vstpat contains harmonics which are frequencies that 
are integrul multiples of the driving frequency. 

In s.me cases the output may centsin cne or more frequencies 
that are 1]wer than the forcing frequency. These frequencies are 
chsracteristically integral submultiples cf the driving function 
and are culled subharmonics. The subharmonics may be small in 
amplitude when compared tu the driving frequency or may be so great 
in comp3riscen thet the forcing frequency may be neglected. 

Unfortunately, there are no general rules available by which 
one may ascertain whether the wccurrence of subharmonic oscillations 
is pussible, Many references un the subject suggest that the question 
Of predicting subharmenics can be answered by considering whether 


oY mvt = simuswidally forced differential equation of the form 


U + k U- 4 ¢éc) = IOOS ypc (2-1) 


LQ 


can have a solution at a frequency vf VU /N, where N is an intezer. 

Equarion (2-1) is a classical equation «fi the theory of non- 
lineezr systems and is knowa as Dutfing’s Equation, 

Tt has been found from experimentation that the occurrence of 
subharmonic oscillations in physical systems is strongly dependent 
on the starting conditions. The amplitude and frequency of the 
driving ferce must fell within certain definite limits, Becsuse of 
this strong dependency on initial conditions it would appear that 
sub-harmonics would only be a transient problem. However, in a 
nunlinesar system the oscillation 18 non-sinusvidal and contains 
harmonics of the fundamental frequency. It is possible to maintain 
the osciilation in a steady-state, under some conditions, by 
supplying esiergy t. the system at any of these harmonic frequencies, 
Hence, the driving force is a harmonic of the fundamental frequency, 
or the fur dsmental frequency is an integral submultiple of the 
fercing frequency. This then is a condition fer subharmonic 


generalion. 


2.1 Duffing’s equation 
Perhaps the most exhaustive study of the subharmonic phenomenon 
ees MO BE on ; 
has been accredited to Haysshi. Ihe following examples of subhar- 
monic oscillation applying Duffing's Equation have been taken from 
his work. 


Corsider the fundsmental equation 
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(2-2) 


in which 246 is a constant damping coefficient ond f(v) is a term 


representing the nenlinear resturing force. 


na 


The period of the forcing function is 2 TY /y, and the subhar- 
monic 1/p has a period 21Y and may be expressed by a linear com- 
bination of sin? and cust . 

What follows 1s an investigation of the relationship between 
the nonlinear characteristic expressed by the term f(v) and the 
urder 1l/v of the subharmonic oscillations. 


Consider the restoring furce to be the polynomial 


wa 2 5 
~ py Co Ge = FO OMCs Cees 1. 
Lr) ile 2. a. (2-3) 
where Cys Cys Cgoeee are constants which are determined by the non- 


linear characteristics, These cunstants are subject to the constraint 


that 
+...= 4 (2-4) 


In stexdy state one may assume that the solution of equation 


(2-2) has the form 


: G Now 
a. rg eR, C + ycos t + EGOS 1 6 





(2227 
‘The approximation for wr is given by 
= B 
Us- = bie 2 (2-6) 


which is legitimate as lung as the nonlinearity is small. 
By substituting (2-5) ints: (2-2) and equating the coefficients 
f the sin and cos? terms separately to zero, the following 
results are obtsined acz.rding to the form of the nonlinear character- 


LS C 1 Co found 1 fe ( 2 > 3) ° 


2 


Case l. f(v) = civ + C3V 


The nonlinearity is symmetrical, or f(v) is an odd function, 
and in this case the constent term z can usually be disregarded, 
Then fer = 254,550.00, @ubstitution as mentioned abeve yields 


(i - 3/4 Cx?4y*)- 3/2 an jy a k a a 


[1-44 (vray) tHeaw'tx - kyo (2-7) 


24 | 
where K= C2 o By simultaneous solution one obtains that 


k(x? + y) = 0 (2-8) 


Eguativn (2-8) implies that the amplitude of subharmonic 
oscillutions is zero as Long as k#0, or damping is present. Hence 
subharmonics of order 1/2, 1/4, 1/5,... camnuet occur in this case. 
But real roots of x and y which do not simultaneously go to zero 
may be obtained for Y = 3, Thereby, one may conclude that sub- 
harmenics cf order 1/3 msy occur when the nonlinear term c,v is 


present. 


Case 2. PCY) =. C2 =e + eae. 
L Z 3 
Since the nunlinesrity is unsymmetrical the constant term z 
in (2-5) must be considered, [Lt is feumnd under rigurous solution 


that subharmonics of the order 1/2 may exist in this case, 


Case 3. f(v) = cay ++ ce 


| ; 3 . ; 
Although the cubic term C3V dwes ick appear in the nonlinearity, 
detsiled investigaticn reveals that subharmonics of order 1/3 still 


occur, along with those of crder 1/5. 


i 


After establishing theoretically the possible occurrence of 
subharmonic solutions to Duffing'’s Equation, Hayashi presents the 


results of several experimental investigstions, a synopsis of which 


follows. 


2el.l Experimental investigation vf subharmonic of order 1/3. 

The experimental data was obtained using an electric oscil- 
latory circult cuntaining a saturable-core inductor and a capacitor 
in series. The circuit equation takes the form of Equation (2-2) 
when a 60 cps a.c. v ltuge is applied tu the circuit. If an initial 
conditiun is prescribed apprepriately, a subharmonic oscillation of 
20 cps msy be started in the circuit. If one theoretically used 


a transformer core whose characteristic is expressed as 


EG re eet ey 


the graph of Figure 2°2 is .btsined. However, the nonlinear 
characteristic of the ordinary transformer core is not truncated 


with the cubic term and is expressed by 


f(y) _ C_V + Cay a Cry. + Cz Vv Ar oa0°0 


ond the region in which 1L/3-harmonic oscillation is sustained 
appesrs in Figure 2el. 

By cormecting a number uf inductance coils in series and ad- 
justing the length of the air gap which is interposed in each core 


% fairly god approximation of the characteristic 
ECV)” “=e + ae 


is ebtauined, 
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Region in which 1/3-harmonic oscillation is 


sustained. 
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Figure 2-1, 
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equation 


When two cores are used, one with an air gap and the other 
without, the experimental verificatii n to the analytical analysis 


is guite satisfacctury. 


2.1.2 Experimental investigation of subhurmunics wef order 1/2. 


Investigsatiuns of systems described by the eqguaticns 





alr cl wom —, “hk ? 
ae fe 7 3 Helos 2a + fs 
At eb (2-9) 
: i ~ 1 - 
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(2-10) 


can be done using an analog computet. Figure 2-3 shows the block 
diegram of an analog computer set up of Equation (2-10). To obtain 
a set up of Equation (2-9) cne needs to replace the servomultiplier 
by an er@inary multiplier which gives a cubic nonlinearity. 

If the initial ccaditions of v(0) and v(0) are chosen appropri- 
ately one may start L/2-harmonic oscillations, Slowly varying the 
values of B and B. produces the regions in which 1/2-harmonics are 
sustained and Figures 2-4 and 2-5 show the experimental results for 
Equations (2-9) and (2-10) respectively. 

Figure 2-6 is a graph of the region in which 1/2-harmonic 


cscillstious are sustained using theoretical calculations. 
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Figure 2-3. Computer block diagram for Equation (2-10). 
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Figure 2-4. Region in which the 1/2-harmonic is sustained; 
nonlinearity by cubic function (analog-computer 
analysis). 
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Figure 2-5. Region in which the 1/2-harmonic oscillation 
is sustained; nonlinearity by symmetrically 
quadratic function (analog-computer analysis). 
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Figure 2-6. Region in which the 1/2-harmonic is sustained 
(calculated). 
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One-half harmonics are also witnessed in the electric circuit 


of Figure 2-7 due to the presence of a saturable-core inductor L 
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Figure 2-/7. Oscillatory circuit containing 
reactor with direct current 
superposed. 


under the impression of an alternating voltage Estné2wt 
The secondary winding is provided on the core so that asymmetry is 
afforded to the nonlinear characteristic by forcing a constant DC 
current flow through it. 

Using the notation of Figure 2-7 we have the equation 


n de + Rep = Escn2wt 


| ; 
Rin = Vig dt 


Leu 


(2-11) 


whete » is the Toumpe Tr of turns in the primary winding and @ is the 


magnetic flux in the core. 


The non dimensional variables u, u_, and v are introduced in 


o 
plece ef i, i_, and @ by che relation 
u 
}2 I ioe I <ee 
i= i rC lo fn A Mh ZS ~ nv | 
(2-12) 
; ae 
where “Ty aad Fr are apprepriate base quantities of the current 


and flux, Neglecting hysteresis the saturation curve has the form 
u -- u a oY sie Cc ay a Ce a Cc va oy ° (2-13) 
The base quantities T, and 9, may be fixed by the relations 


— 3 ’ em é + 5 rp oO = . 
Wy a Cody lees bled C65 4c wae ce 


By eliminating i. and 1 in equations (2-11) one obtains 


ps = 

- a" < e -~ 7 7 = 
de 6 ke dy =p ee Un 9 cae >? © a ~~ © e ¢ pal Rees a AAS & a se 

as Prot 
eft ele 

4 
’ $ ot i ae 
where t = aun = V2 tan K/9. ee fu we 


wee ree Were, 


6 los of + a ae a ac 


In wrder that the experiment more closely egree with the analysis 
a ceompcsite resctwr as described in section 2.1.1 1s used, the 


Satucation characteristics of which are described by 
= a 
(A + Uy s Clu + Ces 


Thus the differentisl equation of the electric circuit becomes 


pe oe 3 “> AS > 
at ot af K de + Ge = ein = fis ce Smmente .. t Oo 
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Figure 2-8 shows the experimental data obtained by varying 


the DC current in the secondary winding and the applied AC voltage. 
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Figure 2-8. Region in which the 1/2-harmonic 
is sustained (experimental). 


ee Subharmonics generated in a closed loop system with thyristors 
amplifiers. 

[It kas been found that the preblem of subharmenics is not 
limited merely to those cases described in the foregoing discussion. 
Dr. F. Fallside and Dr. A. R. Farmer have investigated the eccurrence 
cf subhsrmonics ia systems which employ thyristors to create a 
power amplifier. 

Power amplifiers that use thyristors, thyratrons, magnetic 
amplifiers and mereury-are converters are well known, Like cther 
types of controlled-rectifier amplifiers, the thyristor amplifier is 
a discontinuous element since the input signal to the amplifier con- 
trols the cutput at discrete instants of time only. The cutput wave- 
form 1s characterised by step chianges in voltage and hence is rich 


in cipple. 
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elements preamplifier amplifier elements 
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Fae dire 2-95 General comtrul Leop. 


ae 


When used in a contrel loop such as Figure 2-9, the thyristor 


can be most simply trexzted as a cuntinuvus element, represented by 


a 


its incrementszl DC pain between its steady-state input and its mean 
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output, and the ripple vultage due to the output wave can be dis- 
regarded, This is satisfactory for low bandwidth systems where the 
maximum component frequericy of amplifier input signals is much less 
than the amplifier sampling frequency. Alsw, the large time constants 


incurporated into the lineac filter is sufficient to ensure that 


only a negligible amount of ripple is fed back. 
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Figure 2-10. Waveforms of 2-phase power supply, 
stable operaticn. 


In wide band applications such 2s currest-control loops the 
simple appri<ch is not valid, In such svstems enough g41n must be 
used so that appreciable ripple is fed bsck into the amplifier input 
and, under certszin conditions, 2 self cscillation of the loop, 
known as ripple instability can result. The frequency of the 
oscillation occurs ét en integral submultipie of the basic ripple 
Irequency, i.e., a subharmenic, and manifests itself in the amplifier 
aS a pericdic variation in the firing pattern. 

The wavefurm of Figure 2-10 shows the behavior of the typical 
stable cperation of 4 2-phsse regulsted power supply. 

Figure 2-ll Shows the ripple instability introduced when the 


loop gain is increased by the preamplifier. 
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Pipuwe 2-11. Wavefcurms of 2-phase power supply, presence 
ef cipple irastzbility. 
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Ripple instability cannot be predicted by a simple linearized 
approach to system design and a more accurate analysis must be 


used, 


2.2.1 Amplifier characteristics 
A thytistur smplifier may be generally represented by the 

m=phase atrangement of cells shown in Figure 2-12. The output 
voltage wzveform of the omplifier fer a constant input y is assumed 
t. be as shuwn in Pigure 2-13 fur a resistive wr diode clamped load. 

There ore many types cf firing circuits fer thyristor amplifiers 
and the sinewave reference type is considered here. This type of 
firing has ¢n advantage in that it produces a linear DC gain character- 
istic fur the relation between the amplifier input y and the mean 
output vio leaze ps In this perticular circuit, the input to the 
firing circuit wf each phase consists of the amplifier input plus a 
phased simusdidal reference signal of constant amplitude at the 
excitation frequency. A vate pulse is passed to the thyristor when 
this tutel input passes through zero geing positive. 

The Fallside ard Farmer paper procedes to analyze this parti- 
cular amplifier by two means, an impulse method and a Master 


(scribing Function methiud. 


ee ate nmalvsis by impluse method 
Simce the subharmonic of crder 1/2 is the simpl@st oscillation 
to analyze by the impulse methcd s° that it only will be considered 
here, 
The wibslavced-amplifier cutput waveform characteristic of 


l/2-order subhorm sic instability is shown im Figure 2~l4a. This 
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Figure 2-12. Generalized arrangement of m-phase amplifier 





Bene Wien 24857 Output waveforms of m-phase amplifier, 
resistive or diode-clamped load. 
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a. Unbalanced waveform 


QE 
b. Balanced component 

NE 
c. Pulse train 
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d. Impulse train 


Figure 2-14, Unbalanced operation of a thyristor amplifier 
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waveform is comsidered te have two parts, the balanced component of 
Figure 2-l4b ond the pulse train of Figure 2-14c, alternating at the 
subharmenic frequency of mf. /2. The pulse train which is representa- 
tive 5f the imbalance has no net DC outpuc and injects a signal of 
frequency mf$L/2 into the Lip. A pourtian of this signal, attenuated 
and smocthed by the filter, is fed back into the amplifier input and 
the oscillation msy thus be maintained. Frr eny appreciable oscil- 
lation about the mean firing angle & , the positive and negative 
pulses are different beth in amplitude and duration, Hiwever, for 
diminishingly smaller amounts of imbalance, the pulse train becomes 
in the limit, a symmetrical impulse train of Figure 2-l4c, and 
represents the conditions at the onset of instability. By examining 
the system equations under these conditions ic is possible to formu- 
late a4 gererel conditivn for the loop to sustain a subharmonic of 
1nfinitesimal amplitude. 

Referring to the linesr control loop of Fig. 2-9, assuming 
stable, balanced Operating conditions, the sgmplifier output can be 


written as the Fuurier series 


- S e_(a)sin (nmnet + g) 
ie CCAS n 


(2-14) 
The input tu the chyrister amplifier is 
y = K din) §Ur - Aco | 
(2-15) 
and che equaticn : f the firing circuits at any firing instant is 
iv : ‘ 
y+ Koes ere P)=0 
(2-16) 


Sil 


where 


2irk 
Nt= % + rv) (k20t,2,:- =) (2-17) 
and 
p= Oey r- 
(2-18) 


Substituting for y Equati.n (2-16) beccmes 
W 
K 4) ur, — KY) Alo) + Keos(Nt-F-e)=O —— a-49) 


Considering steady state wperati_n where ve ls constant, the operator 
D may be replaced by jw , after substituting for © , Mt , and iF 
Equation (2-19) becomes 
KA Cjo)u, = K¥GAI4Golu (a) 
< zi’ k 
, ink 
- KD; ze a sin} nim (as mo ) dn + On§ 


n= 


= Yes le =o 


(2-20) 


The operatcr bd) anid Ys, (b) have been replaced by ¥ (;w) and 
aid Y (sw) is the smplituce and On the phase of the expression at a 
frequency W=nmfl, 

If ..e now assumed an infinitesimal amplitude of subharmonic 
oscillation, the system euustivns must be modified by the pulse train 


Av », assuming the pulses have a duration Ax 
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Le = pe oe Sinw sin& cos Im (Nt -% - BLA) 
eee (2-21) 
where 
fs ae (2 Gris a “aie>) 


a2 


| | } Ztr k 
The firing instants now become A= (4+ do) + Vin et iS 


the firing circuit equations bec: mes 


K ¥ Go)u, - K HG AG ao) ee) 


eo ‘ 
Vt , Z my Be Sx WT | 5 d Z iy he 
a K a & ee ee SA mM Sone Cos K m ( “ee 


(=) a 
(2-22) 
shin i ag ge x) +, 8 im S ed Sen} nm(a 
ry 


a we Ot) J +e, -E cos(4-Aa)= O. 


Assuming f\u-eQ , the terms which represent balanced operation 


may be removed by subtracting Eq. (2-20). The conditions by which 


subhsrmonics may exist is seen as 


e ame Aa TT | 
K ; oe Zo om oe SOV rv) Sey A Gee S Qi, 
(: '/2 
cx 3 
2 Z, En(A)nMAK cos(nmas J + ey)f 


An send = O, 


Dividing through by Xda SeNnw 5 and introducing the DC gain 


v& . 
Kaz me eG nny this reduces tx 
TX 


a ZprCn( Xam 
Go P Re CGS o -/ ~ eo ———=~ Cos (nme 
oo hea <4 ' 
O 


(228) 
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This may be written more simply ss 


oO 0} 
Ge} 2% 005% -F2,6,005(4, + Ol} 4 1 =O 
fe im (2-24) 


where GQ) = Sn(alnm 
EX Ky Cos (x-*/m) 


and 
- v= WN A + bn, 


Ta summary, the critical icup gain for ripple instability at 
a subharmenic of order 1/2 con be calculated from the impulse- 
analysis créterion of Eq. (2-24) fur any given thyristor amplifier, 
in any lirear contrel system to any degree of accuracy, by taking 
sufficient terms in the summation. In Eq. (2-24), Cy and C, depend 
the amplifier used and its mean firing angle or cutput voltage; 
the terms Ly oe 


Loop through dC, ww) ¥ (jo) 6 


: O and @. define the remainder of the control 


2.2.3 Describing function for infinitesimal] subharmonic input signals 
cf order 1/2. 

The Fullside and Farmer paper also presents a describing 
function approach tu the analysis of the conditions for sustentation 
of subharmenics of order 1/2, considering the subharmonics appear as 
infinitesimal amplitudes. Again fur a constant input the firing cir- 
Cuif equation is 


ae. Cos Gat ne eee 


(2-16) 
and at the firing instants 
2iVk 
= Cy ee Chan seen 
nC rv , 3 ) (2-17) 


If a sinusbidal signal of infinitesimal amplitude Ay , frequency 
MLL); and some arbitrary phase [A 1S) 2d4easto o. the firing 
equations become 


ry fh 


Vs: M4 cus ( mas cs u) + Costin = Vm) 240 


and ' 
21vkK 
Nt =x tba + EO 


After subtracting Eqn. (2-16), the incremental firing angle 


and the input subharmonic signal amplitude Bu ere related by 


Ay eos 3 9 Coe A gp ZEEE | r Ke 


Zirh fi | 
= = Nex Senos oe ep) | (2-25) 


The addition of a sublsrmonic signal tuo the input of the am- 
plifier in turn priduces «na alternating impulse train at the amplifier 
Guiput. (Fig. 2-14d) This impulse train has a component of sub- 


karmentic frejuency given by 


z2m€ 


mo Bh: 
yp: _ MA ie) 
/2 _ ae a 


ey 
Ten scm 7pm GO. 


Substituting from Eq. (2-25) fur Oy and assuming A& << L 


iy = 2K, cos Yous (GK t +p - #) 


where (2-26) 


Thus, f.r inpset signals of infinitesimal amplitude at subharmonic 
frequeacy a , the amplifier has a gain of 2K, cos and 


intreduces 4 ph2se shift J , 
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The amplifier can then be characterized by the d@scribing 


function, 


{ Vi 
N (5 ") 2 — = Jak cos #f LP 


(2-27) 





cos (Ma we 


4 


This is shown in Fig. 2-15, and is a circle of radius Ka , centered 
| | eof oonry tae 
on ( uo). The normalized version NG =) is a circle of unit 
radius centered on (1, j0). 
This describing function may be used to test for ripple in- 
a" : mr ; 
Stability at the subharmonic frequency = and the procedure can 


be seen in Fig. 2-16. First the open-luop frequency locus KK Do ww}, (yw) 


. 2 . . Ni; Moe 

is plotted; then the normalized describing function N( Mm is 

superimpused, as a unit circle, centered upon the point on the 
mfr : 

frequency locus corresponding to W= “37 . Using conventional 


assumptions, the system will be unstable and oscillate at the sub- 
Ban: | tee ee : Sawee 
hamonic frequency -3= if the circle encloses the point (-1,j0); 


the critical condition being 


tl 
© 


2KKy a 1/2 cos Oa, ele (2-28) 


Referring to Eqn. (2-24), this is the condition which was obtained 


from merely the first term of the exact impulse criterion. 
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Pigubemz—-l>,) Describing function for thyristor 
amplifiers with subharmonic input 
signals of order 1/2. 
Im 






describing function 


Re 


control-loop 
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increasing 


Figure 2-16. Use of describing function to test for 


ripple instability. 
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3. Description of system 

In chapter 2, the SCR amplifier that was discussed was one that 
empluyed a phase firing scheme of control. The thyristors were fired 
individually when the current through each one attained a certain 
phase in relation to a reference sinusoid. The thyristor amplifier 
that is the concern of this paper employs a different method of 
operation. 

The specific system under investigation is a polyphase, thy- 
ristor, AC t> DC rectifier, operating under forced limit cycling 
conditions, The block diagram for the model of the system appears 


in Figure 3= 
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Block diagram fer Regulated Thyristor Power 
Supply. 
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This block diagram my be simplified into the system of Figure 
3-2 where Gy is included to take into account any gain that may be 


assuci-ted with the cumpsarateor, firing circuits and compensation that 
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may be added to the system and is completely linear. G, is the 
linear transfer function of the filter and load, All of the non- 
linegrities of the firing circuit and bridge ore lumped into the 


single describing function N. The nature of this describing function 


e ” ° . i ° 2 e oe 10 
is developed in Leszcynski’s thesis not considering ripple instability. 


a 


Vv + V 
c C o 
A a | 


Premre 3-2, Simplified block diszgram, 


In the phase firing scheme, for polyphase rectifications, it 
is required that there be separate firing circuits for each individual 
thyristor as well as phasing and timing circuits tw ensure proper 
sequence of firing. Hence, 4 three-phase, full wave rectifier would 
require six separate firing circuits, 

So aS te. simplify the cuonutrul philosophy and to curt+il the 
expense vf six sepsrate firing circu1lts, the method considered here 
is t. contre] the entire pulyphase bridge as 4 unit with one firing 
circuit. I: this way the gstes of all the thyristors in the bridge 
are c sneccted in series 21d are pulsed simultaneously with one pulse. 
(A simplified circuit disgrsm can be seen in Figure 3-3.) However, 
only thuse thyristurs whose current has 4a phase between O and 180 


degrees wili fire immedi-.tely os the pulse gces ON. As the phase 


a0 


of the current in the remaining thyristors traverse through 0 degrees 
they too will fire as long as the ON puise is maintained at the 
thyristor gates. Figure 3-4 illustrates the sequence as the six 
thyristors are pulsed ON. As » result, when the firing signal is 

ON, the entire bridge is on and acts as a conventienal diode bridge. 

A slightly different aspect is presented when the bridge is 
pulsed OFF, or the bias is removed from the thyristor gates. When 
the OFF pulse is received, a time delay cccurs before the bridge is 
completely cut vff. This is due to the fact that the thyristor, 
like its gas tube counterpart, the thyratron, is not cut off until 
its shut off counditiuwns have been met, i.e., the bias has been re- 
moved and the collectur vuwltage goes negative with respect to the 
emitter. Figure 3-5 shows the sequence as the six thyristors are 
pulsed OFF. The time delay ic cecurs between the time the OFF 
signal is received snd the last thyristur cuts off, 

A summary of the system operation is as fellows. Assume that 
the desired cutput is a certain voltage and the output voltage is 
beluw this reference level. The rectifier bridge is then supplied 
with an ON sigzal avd the bridge cutput is a three-phased, full-wave, 
rectified voltage which is fed tu the filter. The filter then at-~ 
temusztes the hurm nics “f this veltege and privides a relatively 
smooth DC voltage st the vutput. When the DC output rises above the 
reference level the bridge is pulsed off. Because the shut off con- 
ditions of each thyrist.r must be met, a time delay occurs before 
the bridge turns cff. In the OFF state the input voltage to the 
filiter is esseutislly zero since the reverse current is virtually 


blucked off by a clamping divde. The bridge thus remains in its OFF 
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Figure 3-4. 
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Figure 3-5. Cut-off sequence of thyristor bridge. 


state until the output voltage drops below the reference, whereupon 
the cycle repeats itself. 

Leszcynski’s thesis shows that the output of the power supply 
is a DC voltage with a ripple that is 4 non-sinusoidal limit cycle, 
with a fundamental frequency different from the AC supply voltage or 
any of its harmonics. 

The input tw the filter can be shown as a series of pulses, 


one of which appears in Figure 3-6 belvuw. 


Input to filter (volts) 





ee: | time tx 
Firing pulse y 
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eS | ra st ee t 
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Figure 3-6. Thyristur bridge output pulse 


Suppuse that the period of the supply voltage is T seconds. By 
examining figure 3-6 it can be seen that, depending upon the length 
of the firing pulse, there is a period of T/6, the time that exists 


between one thyristor's firing and the next thyristor’s firing, in 


At 


which the cut off of the firing pulse will not vary the elapsed time 
of the output pulse of the thyristor bridge. It is also seen that 


the time delay c. can vary from T/3 to T/2 seconds in length, or 


Ti3< 6 < T/2 


It is this variation in @ that might give some insight into 
the problem of ripple instability in a forced limit-cycling regulator 


and will be examined more carefully in the next section. 


3.1 Preliminary investigation into ripple instability. 

It was seen in the preceding section that the output of the 
thyristor bridge under steady-state operating conditions is a series 
of pulses that take the form of the pulse shown in figure 3-6. 

If one assumes that the ripple of the DC output of the filter 
is a limit-cycle that closely approximates 4 sinueoid, with a funda- 


mental frequency £, and a fundamental period qT; then the length of 


I 
the firing pulse would be T,/2 or the length of time that ie is less 
than es in one limit-cycle period T- 

Figure 3-7a is a replica of figure 3-6 in which the firing 
pulse is T,/2 seconds Long anda the thyristor bridge pulse length is 
ty seconds. Figure 3-/7b shows the firing pulse still to be T,/2 
seconds but displaced in time from the first firing pulse by€ 
seconds. The resulting pulse formed by the thyristor bridge is t, 
seconds which can be seen to be less than tie 

Suppose one used this sine wave limit cycle over a long period 
of time to generate a train of puises that might be obtained from 


the thyristor bridge under steady-state operating conditions. Assume 


that the bridge is pulsed on when the sine wave is greater than its 
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Figure 3-/b. 


46 


‘asyTng JO AsqunNNY SA spuodsas UT YIPTM es {Nd °9-¢ dINSTY 
°Zq19y Ty = aTOAD-FIUTT 8° q 
0Z OT 0 











Tone a 
a 





Tz20°0 


47 





ft 


‘as[ng FO 1aqUMN SA spuodas UT YIPTM as{ng °g-¢ ain3Ty 


°z3zey LZ = aToA0-ITuUTYI = *p 


—TARABANIAL 
a 


"23194 €€ = VTIAD-FIWTT =o 


MAA ll 
Ml bilialiee > — 2 





£c0°0 











48 


asTNg JO AaqUMN SA Spuodag UT YAPTM es [Ng °9-¢ o1n8Ty 


*Z9294 GZ = VTOAD-JIWTI a) 





49 


ue and is shut vff when the sine wave gues below its mean, 
sc as to simulate the freqpency of the errcr signal becoming greater 
or less than zero. Observed cver fifty or sixty cycles of the 


frequency f, one notes that the pulse train from the bridge takes 


1 
wen an obvious pattern as cam be seen in figure 3-84-e. 

Observe that the pericd of repetition of the pattern established 
by this means is some times many times greater than the Limit cycle 
period T, but im all cases is an integral multiple of the limit-cycle 
period, If vwne were tu write « Fuurier Series representation of 
this pulse train, as complicated as it might be, the fundamental fre- 
quency in each csse would be a subharmonic of the sine wave that 
generated it. Hence, it can be seen that a pulse train generated by 


@ simusvidal signal with frequency f, presents energy to the filter 


1 
that contains subharmonics of the generating signal. 

This, of ccurse, is not irrevocable proof that subharmonics do 
exist in the output because it is known that the Llimit-cycle is not 
da true sinusvid and also one does not know what form of pulse train 
is generated by a2 Limit-cycle in which ripple instability does exist, 
but one must admit to a definite possibility that subharmonics of 
the Limit-cycle frequency might be present in the cutput. 

The mexc step would be to simulate the system on an analog or 
digitsl computer, if nut the actual plant itself, and in fact see 
if subharm nics sre present in the output and under what conditions 


they exist. This is indeed what will be dune in the next chapter. 
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4. Synthesis of the forced limit-cycling regulator. 

The system under study was successfully synthesized in 
Leszcynski's thesis by a digital computer program utilizing Fortran 
IV language and the Naval Postgraduate School IBM 360 Computer Facility. 
In the program the firing circuit and SCR bridge was simulated by 
logic statement manipulation and the output of the linear filter was 
obtained by calling a library subroutine INTEG 2 which gives a Runge- 
Kutta solution to a group of differential equations, in this case 
the state variable equaticns developed from the dynamics of the linear 
filter (see Appendix 1). The linear filter is assumed to be as shown 
in Figure 4-1. The transfer function for the filter and load is 


developed as follows: 





Ua = Cie a C= aig, + oe 
‘ din 
= _ 4+ Ue, Te ss Ve 
7 dt At RR 
ace ! ( ~ %) Aux ee ~ ies 
ee a : d+ @ [eres 


At 


whe re Uo = L% 


Choosing the state variables as 
K, = UG x Ya fe 7 K. {re 


Keay a Ceo rg K/L 


a signal flow graph of the linear system is shown in Figure 4-2. 


Using Mason's gain rule the transfer function in Laplace 


ee 


7 
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NOLALLON is: 


at 


a, 


in 


Figure 4-1. 


Figure 4-2, 


Filter and load, 


Signal flow graph for linear system. 
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out 





The feedback path between the error signal und the firing 
circuit is known in the actual system to contain a dead zone. Figure 
4-3 1s the dead zone that was assumed in simulation of this system, 
with O& being the width of the dead zune. In the program a gain C(6) 
is intreduced into the feedback path tu provice amplification to the 
error signal prior to the dead z.»ne and has the ability of defeating 
the effects of the dead zune if used properly as will be seen in the 
experimental results. 

In conducting the experiment the supply voltage was considered 
to be a three-phase 60 cycle AC source with a peak amplitude of 65 
volts. The value of the inductance and capacitance of the Linear 
filter were kept constant throughout the experiment at 0.0005 henries 
and 0.1 farads respectively. 

Thus, the experimental data was collected using the following 


FaCLOrS as variables: 


C(4) - The value of the load resistance in ohms 


C(5) - The value of the reference voltage 


i 


C(6) - Ihe amount of error signal amplification 


C(7) - The width of the dead zone (A ) 


Q 


Using a step size of 9.00002222 secunds it was found that in most 
cases by setting the initial conditions clese to the desired reference 
voltage the transient died cut within the first 0.1 second and the 
effect of ripple instability could be successfully observed from 0.1 
tu 0.2 seconds. Thus, the graphs that are shown begin at 0.097687 
seconds and continue to 0.19702 seconds. The effects of manipulating 
the several variables can well be seen during this time frame in 


most cases, 


520, 


Signal out Signal out 





NOTE: The shaded areas 
denote the period in which 


the SCR bridge is ON, 


Signal in 


Figure 4-3. The dead zone characteristics of the feedback 
system. 
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4.1 Experimental results 

The first concern was te find wut if the output ef the filter 
actually contained some subharmonic ripple, so taking time intervals 
of 0.1 seconds the simulation program was run until the transient 
had obviously died sway and the steady state cutput was examined, 
Fig. 4-4b. shows the steady state cutput of the regulator with the 
load equal t. 1.6 chms and the reference voltage ser ar 60.0 volts. 
(Ihe desd zone in esch case is set at 0.0 vwlts unless ctherwise 
nected on the graphs.) 

The load resistance was subsequently doubled and then halved to 
see the respective resuits on the output. Those results can be seen 
in Fig, 4-4a and c,. Fig. 4-5 shcws the input pulse train for the 
correspending time period. 

In Fig. 4-62 the reference voltage with the load set at 3.2 
chms was drcpped to 58.0 volts. Fiz. 4-6b presents the results 
when the errer signal was doubled but zeru dead zone was considered. 

Keeping the load resistance at 0.8 ohms, the reference voltage 
was dripped in steps from 65.0 veits to 48.0 volts. These results 
appear in Figs. 4-8, 4-9, and 4-10. 

Using the values of R=0.8 hms and ae =60.0 volts the system 
wa$ exemined assuming the dead zcne tu be 0.5 velts. These results 
were compsred tv the results cbtained by setting V of 00«? volts and 
the dead game set at 0.0, Fig. 4-11 shows this comparison. 

Figure 4-13 is a comparison between the results ubtained by 
setting Vig f 02-0 volts and the dead zune at 0.0 and setting Vif 0009 
velts and dead zone equal te 2.0 volts. The errér signal was then 


amplified by 2 factor wf 1.5 and compared to the praceeding experiments. 
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Figure 4-12. Vz. 
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Finally the system was examined when the pulse train was 
amplified by a factor of 2.0 befwre going into the filter and the 
results were compared to setting the source voltage to double its 


original value. These experimental results appear in Fig. 4-14. 


4,1 Summary and analysis of the experimental data. 
The transfer function for the linear filter as stated previously 
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In the Bode furm this becomes, 
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In decibels, 
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Comparing this to the n-rmel form cf the Bede plut equation 
T(jW) = +20 Log | (jM/y 4 oan (j &) ) ae | 
db n n 
it can be seen that 
= ee 28 fy, =U/R and © = 1/2R 
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Using these three values of load resistance the Bode diagrams 


for the linear filter response con be seen in Figure 4-16. 


Table 4-1 is a tabular summary of the results of the experiment 


in which the reference voltage was kept censtant at 60.0 volts and 


the value of the load resistance was varied. In Table 4-1, 
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Frequency of the fundamental lLimit-cycle in 
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Frequency of the ebvious subharmonic ia 
hertz 


Peak to peak amplitude of the larger spike 
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Peak t peak amplitude of the smallest spike 
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Damping coefficient 
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uss 5.904 4.192 044 


were ae ere gp PS mere me ee 8 ee 6 ee res gt TT ree ca | 
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The values of the subharmonic frequencies are noted for each 
case on the Bode diagrams of Figure 4-16, 

Table 4-2 is 3a tabulsted summary of the experiment in which the 
load resistence was held cunstant st 0.8 ohms and the reference volt- 


age was decreased. 





Ror 0g8eOhms 





Veet one *SuB_ i 
65.0 27.56 8 10.004 
64.0 32.91 0 8.001 
62.0 35.79 17.89 6.927 
60.0 45,06 22 5.904 
58.0 51.79 25.89 5.003 
56.0 60.13 15.03" 4.997 
54.0 51.79 36. 36 32301 
52eoe 60.93 30.47 4.410 
50.0 90.25 45.13 2.592 
99 .99 2 ere 


af, 


wate val 


| 
| 48.0 


hpie 4-2. 


1/4-order subharmonic 


" Data does not fullew pattern. 





Because the physical system is known tu contain 4 dead zone in 
the feedback path, Table 4-3 is 2 compariscn of the system's response 
when a dead zone is assumed and the system's response when no dead 


zone is cwnsidered, 
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Figure 4-16. Response of Linear Filter 


Gain in(db's) vs. frequency in radians/sec. 
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Ji Conclusicns 

The main conclusion that cen be drawn by examination of the 
experimental data is that under certain conditions the forced limit- 
cycling regulator discussed will support subharmonic ripple instability. 
The other object of the investigation, that of determining the exact 
causes for the subharmonic sustentation, was not clearly achieved 
and the cunclusions that might be drawn toward this end are somewhat 
Speculative and will appear in the next section in a discussion of 


the recommendations for further study in this area, 


6. Recommended areas fur further study 

One tentative conclusion that might possibly be reached in 
reviewing the experimental results is that the subharmonic that is 
sustained in the system is dependent upon the numerics of the com- 
ponents of the filter that is employed. The filter in this case was 
both simple and a little unrealistic but it did bring the interesting 
observation that the subharmonic that was the most evident in the out- 
put was one that was quite near the resonsnt frequency of the linear 
filter, Examining the results presented in Figure 4-4 one can see 
that althowgh the damping of the linear filter is increased as the 
load resistance is Lowered, the amplitude of the subharmonic is in- 
creased, This might be explained by the Bode diagram of Figure 4-16 
where it can be seen that as the load is reduced the limit-cycle 
frequency is also reduced and its 1/2-order subharmonic is moved 
closer to the resonant frequency of the filter. 

In Figure 4-9b, where the resistance is 0.8 ohms and the 
reference voltage is 56.0 velts, the 1/4-urder subharmonic of the 
limit-cycle is more evident than the 1/2-order subharmonic as was 
the case with the rest of the results. This also could be due to 
the fact that this crder subharmenic was closer to the resonant peak 
than was the 1/2-.rder subharmonic. 

Based on the above philosophy, a recommendation for further 
study would be to refine and design the filter to produce a desired 
Limit-cycle whose subharmonics are well attenuated by the Linear 
Paitere 

it can be noted in Table 4-3 that as the dead zone in the 


feedback path was increszsed the limit-cycle frequency decreased and 
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its 1/2-order subharmonic mcved eway from the filter's resonant peak. 
Any further work should also include a mre detailed study of the 
effects of the dead zone on the system. 

Another area fur study that was touched on briefly but in no 
extensive detail was the matter of varying the supply voltage and 
studying its effect on the ripple instabiliry. 

In the realm of mathematical analysis several of the references 
noted in the bibliography suggest that a possible method of predicting 
the perfurmance of the system would be to employ a Dual Input 
Describing Function or even a Multiple Input Describing Function. 
Surely, any further study should give some consideration to these 
analytical methods in attempting to understand ripple instability in 


this system. 
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APPENDIX I 


Computer program for simulation of system with numerical data 


for the conditions when the reference voltage is held constant at 
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Some Nonlinear Feedback Control Systems under high gain conditions exhibit 
the phenomenon of subharmonic instability, or cuntain subharmonics of the 
fundamental output frequency. A general discussion of subharmonics in nonlinear 
systems is followed by an investigation of ripple instability in a forced limit - 
cycling voltage regulator containing a thyristor or SCR bridge utilizing an ON- 
OFF switching scheme. 


A digital computer program is used to simulate the dynamic response of 
the system under different loading conditions and for different reference voltage 
levels. 


ID Fok".1473 (Pace 1) UNCLASSIFIED 


iN 0101 -807-6811 Security Classification A~31408 


UNCLASSIFIED 


~ Security Classification 


14 
NI A 











Ripple 
Limit~cycle 
Subharmonics 


Thyristor 





DD er al 4 73 (BACK ) UNCLASSIFIED 


S/N 0101-807-6821 Security Classification a> Sa 





162 























